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Abstract 

Using a large sample of charm semileptonic decays collected by the FOCUS pho- 

toproduction experiment at Fermilab, we present evidence for a small, even spin 

K~7r'^ amplitude that interferes with the dominant K component in the — > 

K^TT^ fi^iy final state. Although this interference significantly distorts the D"*" — ^ 

K^TT^ decay angular distributions, the new amplitude creates only a very small 

distortion to the observed kaon pion mass distribution when integrated over the 

— *o 

other kinematic variables describing the decay. Our data can be described by K 
interference with either a constant amplitude or broad spin zero resonance. 
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This paper describes discrepancies between the observed decay intensity for 
the decay K^n^n~^u and that expected for pure K yU+z/ decay. 

The data are overlayed with a simple model incorporating a particular choice 
of interference amplitude. A later paper will present the results of fits for the 
parameters describing the interfering amplitude along with new values of the 
form factors describing — > K jji^v. 

Five kinematic variables that uniquely describe K~tt^ ii^v decay are 

illustrated in Figure 1. These are the K^tx^ invariant mass {mKir) , the square 
of the fjLu mass (q^), and three decay angles: the angle between the tt and the 
D direction in the K~7r~^ rest frame (^v), the angle between the u and the 
D direction in the yuz/ rest frame (Oe), and the acoplanarity angle between the 
two decay planes (%). The acoplanarity conventions for the and D~ will 
be discussed near the end of this paper. 

It has been known for many years that the L>+ —>■ K~t:^ ii^u final state is 
strongly dominated by the K fi'^u channel [1]. In the process of study- 

ing the K^Ti^fi^u decay distribution, we found significant discrepan- 

cies between our data and the angular decay distributions for — > K fi'^v 
predicted using previously measured [2] form factor ratios. In particular, we 
discovered a significant forward-backward asymmetry in suggesting a new 
term in the decay distribution that is linear in cos 9\ and a strong function of 
the kaon-pion mass. In pure K /i"*"//, the acoplanarity-averaged decay 

intensity should consist of terms with only even powers of cos In this paper 
we present an explicit model for this effect in terms of coherent interference 
between a small spin zero component of the kaon-pion system with a dominant 
— > K n'^u component. Throughout this paper, we describe this interfer- 
ence as due to a constant, s-wave amplitude described by a single modulus 
and phase. We do not exclude the possibility of a broad spin zero resonance 
or even higher spin resonances with helicity amplitudes tuned to resemble the 
behavior of a single, broad s-wave resonance. 

Angular momentum considerations allow us to predict the dependence of the 
interference-induced, linear cos^^v term on other decay variables such as g^, 
9g , and %. We present this model and show how well our data match these 
predictions. Throughout this paper, unless explicitly stated otherwise, the 
charge conjugate is also implied when a decay mode of a specific charge is 
stated. 

The data for this paper were collected in the Wideband photoproduction ex- 
periment FOCUS during the Fermilab 1996-1997 fixed-target run. In FOCUS, 
a forward multi-particle spectrometer is used to measure the interactions of 
high energy photons on a segmented BeO target. The FOCUS detector is a 
large aperture, fixed-target spectrometer with excellent vertexing and particle 
identification. Most of the FOCUS experiment and analysis techniques have 
been described previously [3]. Here, for the first time, we describe our muon 
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Fig. 1. Definition of kinematic variables. 



identification. FOCUS uses two muon detector systems, the inner muon ho- 
doscope, and the outer muon detector. Muons are identified by their abihty 
to penetrate approximately 21 interaction lengths of absorber for the inner 
muon system, and 18 interaction lengths for the outer system. For the inner 
system, potential muon tracks are projected through the electromagnetic and 
hadronic calorimeters, and additional iron shielding walls. This trajectory is 
then matched to hits recorded in an inner muon detector consisting of six ar- 
rays of scintillation counters subtending approximately ±45 mrad. The outer 
muon detector consists of three views of resistive plate chambers which are 
shielded by the outer electromagnetic calorimeter and the iron yoke of the 
second analysis magnet. This subtends an additional region of roughly ±140 
mrad. For the outer muon system, potential muon tracks are projected through 
the magnet yoke and then matched to outer muon hits. 

Our analysis cuts were chosen to give reasonably uniform acceptance over the 5 
kinematic decay variables, while still maintaining a reasonably strong rejection 
of backgrounds. To isolate the K^-n-^fi^u topology, we required that 

candidate muon, pion, and kaon tracks appeared in a secondary vertex with 
a confidence level exceeding 5%. The muon track, when extrapolated to the 
shielded muon arrays, was required to match muon hits with a confidence level 
exceeding 5%. The kaon was required to have a Cerenkov light pattern more 
consistent with that for a kaon than that for a pion by 2 units of log likelihood, 
while the pion track was required to have a light pattern favoring the pion 
hypothesis over that for the kaon by 2 units [4]. 

To further reduce muon misidentification, an inner muon candidate was al- 
lowed to have at most one missing hit in the 6 planes comprising our inner 
muon system. In order to suppress muons from pions and kaons decaying in 
our spectrometer, we required inner muon candidates to have an energy ex- 
ceeding 8 GeV. For outer muons we required an energy exceeding 6 GeV. 
Non-charm and random combinatoric backgrounds were reduced by requir- 
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ing both a detachment between the vertex containing the K~'k'^ii^ and the 
primary production vertex of 12 standard deviations and a minimum visible 
energy {Ek + E^, + E^) of 30 GeV. To suppress possible backgrounds from 
higher multiplicity charm decay, we isolate the Kt:^ vertex from other tracks 
in the event (not including tracks in the primary vertex) by requiring that 
the maximum confidence level for another track to form a vertex with the 
candidate be less than 0.1%. 

In order to allow for the missing energy of the neutrino in this semileptonic 
D'^ decay, we required the reconstructed Kttii mass be less than the nominal 
mass. Background from — > X^tt+tt'*", where a pion is misidentified as a 
muon, was reduced using a mass cut: we required that when these three tracks 
were reconstructed as a Kinr, their Knir invariant mass differed from the 
nominal mass by at least three standard deviations. In order to suppress 
background from D*~^ —>■ D^-n^ —>■ {K^ u)ti^ we required M{K^ n^un^) — 
M{K~ ijf^v) > 0.18 GeV/c^. The momentum of the undetected neutrino was 
estimated from the line-of-flight as discussed below. 

We assumed that the reconstructed D momentum vector points along the 
direction defined by the primary and secondary vertices. This leaves a two- 
fold ambiguity on the neutrino momentum. We use the solution that gives the 
lower D momentum, which, according to our Monte Carlo studies, produced 
somewhat better estimates for the kinematic variables. Due to resolution, 50% 
of events were reconstructed outside physical limits (the pj_ of the charged 
daughters relative to the direction implied a parent mass larger than the 
nominal mass). These events are recovered by moving the primary vertex 
to the nearest allowed solution and the kinematic variables are recomputed. 
Monte Carlo studies show that the inclusion of the recovered events does not 
significantly degrade either the resolution or the signal-to-noise ratio. 

It was important to test the fidelity of the simulation with respect to reproduc- 
ing the resolution of those kinematic variables which depend on the neutrino 
momentum. To do this, we studied fully-reconstructed D° — > K'n'^Tr'^n'' de- 
cays where, as a test, one of the pions was reconstructed using our line-of-flight 
technique. We then compared its reconstructed momentum to its original, 
magnetic reconstruction in order to obtain an "observed" resolution function 
that was well matched by our simulation. 

Figure 2 shows the K~7r'^ mass distribution for the signal we obtained using 
the cuts described above. A very strong K (896) signal is present. To as- 
sess the level of non-charm backgrounds, we plot the "right-sign" (where the 
kaon and muon have the opposite charge) and "wrong-sign" K~7r^ mass dis- 
tributions separately. We subtract the distributions of wrong-sign from right- 
sign events as a means of subtracting non-charm backgrounds that are nearly 
charge symmetric. 
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Fig. 2. D"*" ^ K~TT~^ fi'^L' signal. Right-sign and wrong-sign samples are shown. The 
approximate wrong-sign-subtracted yield is 31 254 events. In the mass window from 
0.8-1.0 GeV/c^ there is a right-sign excess of 27 178 events. A Monte Carlo that 
simulates the production and decay of all known charm species predicts that ^7% 
of this excess is actually background from other charm decays. 

Figure 3 compares the wrong-sign-subtracted cos^^v distribution to that pre- 
dicted by our Monte Carlo which incorporates all acceptance, and resolution 
effects as well as all known charm backgrounds. We show separate distribu- 
tions for events with a K~'ir~^ mass in the range 0.8 < rriKn < 0.9 GeV/c^ and 
0.9 < niKn < 1.0GeV/c2. 

Although the events with a reconstructed mass above 0.9 GeV/c^ are a rea- 
sonable match to the prediction, a striking discrepancy is apparent in the 
cos6'v distribution for those events below the pole. As will be explained in 
detail later, we believe this can be explained by the interference of a broad 
(or nearly constant) s-wave amplitude with the Breit-Wigner amplitude de- 
scribing the K . In particular, this interference creates a term in the decay 
intensity that, when averaged over acoplanarity, is linear in cos^v whereas 
the pure K jX^v process produces only even powers in cos^v in this 

intensity. The slight forward-backward asymmetry present in the Monte Carlo 
histograms of Figure 3 primarily reflects acceptance variation.^ 

We exploit the fact that all acoplanarity-averaged terms in the decay intensity 
expected for pure — > K /i+i/ should be proportional to even powers of 
cos6'v (see Eqn. 1 below). Because our analysis cuts give reasonably uniform 
acceptance over cos6'v, we can construct a weight designed to project out any 
Unear cos 6'v contribution to the decay distribution. This weight is the product 



^ The known charm backgrounds tend to have the opposite cos 6\ asymmetry to 
that observed in Figure 3. Their level can also be significantly reduced under tighter 
analysis cuts that still preserve the much larger asymmetry in the data. 
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Fig. 3. Event distribution in cos^V) split between samples above and below 0.9 
GeV/cP. The points with error bars are (wrong-sign subtracted) FOCUS data and 
the solid histogram is a Monte Carlo simulation, including the signal with the mea- 
sured form factor ratios [2] and all known charm backgrounds. The Monte Carlo is 
normalized by area for each plot independently. 

of a wrong-sign subtraction weight (+1 for right-sign and -1 for wrong-sign) 
multiplied by cos6'v. The weighted mK-K distribution is shown in Figure 4 and 
compared to two Monte Carlo simulations. One Monte Carlo is based on pure 
— > K /i+z/ and known charm backgrounds (dashed histogram) while the 
second (solid histogram) also includes the interfering s-wave amplitude that 
we will describe later. The dashed histogram shows that the residual effects 
of charm backgrounds, acceptance variation and resolution effects produce a 
much smaller variation with cos 6y than we observe. The striking mass depen- 
dence of the linear cos dy term displayed in Figure 4 will be the principal tool 
we will use to estimate the parameters of the interfering amplitude. Figures 
5-6 compare the dependence of the cos 6y asymmetry on two other kinematic 
variables (cos 6'^ and g^) to the Monte Carlo with and without the s-wave am- 
plitude. The acoplanarity dependence of the interference term will be discussed 
in a later section. 

It was possible to understand the forward-backward asymmetry in cos 9y using 
a simple model summarized by Eqn. 1. Using the notation of [5], we write the 
decay distribution for K^-K^fi^u in terms of the three helicity basis 

form factors: if+ , Hq , For simplicity, we show the decay distribution in 
the limit of zero lepton mass 



drfiK-K dq^ dcosOy dcosOg dx 



oc Kq'^ 



(1 + cos^^) smeve'^BK*oH+ 

- (1 - cos^^) sinOye-'^BK^oH^ 

- 2 sin ^^(cos 9yBK*o + Ae'^)Ho 



where K is the momentum of the K tt^ system in the rest frame of the D~^. 
The helicity basis form factors are given by: 
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(Mfl + mx.)Ai(g2)^2 
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{Ml - ml 



g2)(Mz5 + m^,)A(g2) -4 



MlK^ 




The vector and axial form factors are parameterized by a pole dominance 
form: 



where we use world average [2] values of Ry = V{0)/Ai{0) = 1.82 and R2 = 
A2(0)/Ai(0) = 0.78 and nominal (spectroscopic) pole masses of Ma = 2.5 GeV/c^ 
and My = 2.1 GeV/c^.Q 

The B^'.o stands for a Breit-Wigner amplitude (Eqn. 2) describing the K 
resonance: 



In Eqn. 1, the s-wave amplitude is modeled as a constant (no variation with 
mKn) with modulus A and phase S. Angular momentum conservation restricts 
its contribution to the Hq piece that describes the amplitude for having the 
virtual in a zero helicity state relative to its momentum vector in the 
rest-frame. 

Assuming this new, previously unreported amplitude is small, it will primarily 
affect the decay distribution through interference with the dominant Breit- 
Wigner amplitude. Expanding Eqn. 1, we find that interference between the 
s-wave amplitude and the B^-o amplitude produces the following three inter- 
ference terms: 



^ Eqn. 1 implicitly assumes that the q'^ dependence of the s-wave amplitude coupling 
to the virtual W'^ is the same as the Hq form factor describing the — > K*^ fi^u. 
This assumption is consistent with our data as illustrated in Figure 6. The modulus 
A would then be the form factor ratio. 

^ We are using a p-wave Breit-Wigner form with a width proportional to the cube 
of the kaon momentum in the kaon-pion rest frame (P*) over the value of this 
momentum when the kaon-pion mass equals the resonant mass (-Pq). 
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Interfere = 8 cos sin^ OeA^ (e-^^B^^.o) (3-a) 
- 4(1 + cos Oe) sin 6^ sin (e^^^-'^^B^-o) H+Ho (3-b) 

+ 4(1 - cos Oi) sin Oe sin (e^^^^+^^B^^.o) H^Hq (3-c) 



Hq, will be present if 



Only the first term (3-a), 8 cos sin^ ^^^3? (e-^-'Bj^^o 
one integrates over the acoplanarity variable x- Since our acceptance is very 
uniform in x, we will primarily observe the effects of 3-a when we average 
over X- We will begin by studying the acoplanarity-averaged asymmetry dis- 
tributions before turning attention to the two acoplanarity dependent terms: 
Eqn. 3-b and 3-c. 
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Fig. 4. Asymmetry distribution in Kit invariant mass, (a) The dashed line represents 
our Monte Carlo simulation with no interfering s-wave amplitude. The experimental 
data are the points with error bars. The solid line is the Monte Carlo with an 
s-wave amplitude of approximately 0.36 (GeV)^^, and a phase of j. Known charm 

backgrounds are simulated for both, (b) A plot of 3ft (^e~^^Bj^*o^ versus itiktv for three 
choices of the phase 6. We also show an alternative modeling of the s-wave amplitude 
as a broad (F = 0.4 GeV/cP) resonance with a mass of 1.1 GeV/c^. We have put the 
broad, s-wave resonance in with a real phase relative to the K*^ Breit-Wigner, as one 
might expect given the presumed absence of final state interactions in semileptonic 
decay. This resonance solution is not unique. 



The previously mentioned, weighted mKn distribution is shown in Figure 4. 
The shape of the cos term versus m^Tr is a strong function of the interfering 
s-wave amplitude phase 6. Because acceptance and resolution corrections are 
small, the phase can be informally determined from the rriKn dependence 
expected from the interference of this phase with the phase variation expected 
for a Breit-Wigner, *o ) . Figure 4(a) demonstrates that the cos 6y 

weighted distribution in data is consistent with a constant s-wave amplitude 
of the form 0.36 exp(z7r/4) (GeV)~^. The magnitude of this amplitude is 
the value required to match the total asymmetry in data over the interval 
0.8 < mK-K < 0.9 GeV/c^. In the discussion to follow, the s-wave amplitude 
will be fixed to this value. 
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We next turn to a discussion of the dependence of the acoplanarity-averaged 
hnear cos^v term on other kinematic variables. Figure 5 compares the ob- 
served cos 6'v weighted distribution as a function of cos 6'^ to that expected in 
our model using our constant s-wave amplitude of 0.36 exp(z7r/4) (GeV)~^ 
We show the distributions for m^Tr both above and below 0.9 GeV/c^. 

In the context of our model, the acoplanarity-averaged, linear cos^v term 
should be of the form: 8 cos^v sin^ ^^^43? (^e~*'^Bj^.o^ ifp and is proportional 
to 1 — cos^ 6 1. This parabolic dependence is quite evident in the weighted 
histogram shown in Figure 5 for those events with mK-K < 0.9 GeV/c^ where 
the forward-backward asymmetry is the largest. 
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Fig. 5. Asymmetry distribution in cosOi below (a) and above (b) the K*^ pole. The 
dashed line represents our Monte Carlo simulation with no interfering s-wave ampli- 
tude. The solid line is the Monte Carlo with an s-wave amplitude of approximately 
0.36 (GeV)~^, and a phase of ^. We expect the acoplanarity-averaged interference 
to be proportional to 1 — cos^ and to appear predominantly below the pole. Be- 
cause the cos coefficient is negative, the weighted distribution appears inverted. 
The model with the s-wave amplitude is in good agreement with the data. 

As a final test of the acoplanarity-averaged interference term, we examine 
the dependence of the linear cos^v coefficient. In our model, this coefficient 
should be proportional to K Hl{q^). Figure 6 compares the cos 6'v weighted 
q^ distribution to the data with and without the additional s-wave amplitude 
of 0.36 exp(i7r/4) (GeV)-^ 

Eqn. 3-b and 3-c produce an s-wave interference term with sinusoidal variation 
in X s-iid are proportional to sin 6'v. In the absence of the s-wave amplitude, 
all acoplanarity dependences either involve cos(2x) or are odd functions of 
cos^v The fact that (with the known form factors) the H_ term (Eqn. 3-c) 
dominates over the term has guided us in the construction of a weight to 
study the acoplanarity dependence of the s-wave interference. We weighted the 
acoplanarity distributions by weights of the form sin(x+(5) and cos(x+5) where 
(5 = 7r/4 (the phase of s-wave amplitude in our model). We also multiplied 
these weights by for right-sign events and —1 for wrong-sign events to 
subtract away the bulk of the non-charm background. For reasonably uniform 
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0.0 0.2 0.4 0.6 0.8 1.0 
GeVVc* 

Fig. 6. Asymmetry as a function of . The solid curve is our Monte Carlo in- 
cluding an interfering s-wave amplitude of 0.36 exp(i7r/4) (GeV)^^. The dashed 
curve has no interfering s-wave amplitude. Both Monte Carlo simulations include 
known charm backgrounds. The rightmost bin also includes those few events where 
> l.l(GeV/c2)2. The model with the s-wave amplitude is in good agreement 
with the data. 

acceptance, these weights should average to zero for any constant or cos(2x) 
terms or any of the cos^v cosx terms present without s-wave interference. 
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Because the H_ term (Eqn. 3-c) dominates over the if+ term, by offsetting the 
phase of the weight by the phase of our s-wave amphtude, we have arranged 
things such that the cos(x -|- 6) weighted distribution should be proportional 
to 3?(B/^*o) (an odd function of 171x^,-1710) and the sin(x-|-5) weighted distribu- 
tion should be proportional to $j(B^.o) (an even function of rriK-K-mQ). These 
expectations are essentially borne out in the weighted plots shown in Figure 7. 
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Fig. 7. Test of acoplanarity interference terms. Plot (a) uses the cos(x + <5) and is 
expected to be proportional to the real part of a B^*o amplitude; while the plot (b) 
uses the sin(x + S) weighting and is expected to be proportional to the imaginary 
part of a B^»o amplitude. In both figures, the solid histogram is our Monte Carlo 
including the 0.36 exp(i7r/4) (GeV)~^ term; while in the dashed histogram there is 
no s-wave amplitude. 
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The presence of s-wave interference creates an interesting new complication 
in the conventions concerning the definition of the acoplanarity variable x- 
Our conventions are guided by CP symmetry considerations. We define the 
sense of the acoplanarity variable via a cross product expression of the form: 
(P^ X Py) X [Pk X Ptt) ■ Pktt where all momenta vectors are in the rest 
frame. Since our x convention involves five momenta vectors, we believe that 
as one goes from D~ one must change x ~^ ~X Eqn. 1. In the 

absence of the interference, there is no need to consider the sign conventions 
on X since the decay distribution involves only cosines of x or 2x- 



This point is made graphically by Figure 8 which compares the observed 
wrong-sign subtracted acoplanarity distribution for the and D~ . These 
distributions are quite consistent once the sign convention is properly reversed 
for the D~ relative to the , indicating no evidence for CP violation in these 
decays. 
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Fig. 8. The wrong-sign subtracted acoplanarity distribution separated by charm. 
The "x" points are for the while the "diamond" points are for the D~ . (a) 
compares the distributions without the required change in the x convention as 
discussed above, (b) compares the distributions with the correct x sign convention 
change. 



We have presented compelling evidence for the existence of a coherent K^ti^ 
s-wave contribution to K^n^ ^^u. It has been assumed in all previous 

experimental analyses that this decay was strongly dominated by the process 
K yU"*"//. The previously unobserved, s-wave contribution is modeled as 
a constant amplitude of the approximate value 0.36 exp(i7r/4) (GeV)~^. Its 
strength, 0.36, is roughly 7% of the K Breit-Wigner amplitude at the pole 
mass in the term that couples to Hq in Eqn. 1. The effect of this new inter- 
ference is very noticeable in our data and creates a ~15% forward-backward 
asymmetry in the variable cos^v for K jji^v events with rriK-K below 

the K* pole. 



Although such an interference effect has been discussed in the phenomenolog- 
ical literature [6] , there has been no discussion of it (to our knowledge) in the 
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experimental literature. How could such a large effect have gone unnoticed in 
the past? We believe one answer is that an amplitude of this strength and form 
creates a very minor modulation to the mK-K mass spectrum as shown in Fig- 
ure 9. Another reason is that this effect is much more evident when one divides 
the data above and below the K pole as we have done. We were unable to find 
evidence that this particular split was studied in previously published data. 
Finally, the FOCUS data set has significantly more clean K~tt^ 
events than previously published data. 




Fig. 9. We show the wrong-sign subtracted niK-w distribution in data (points with 
error bars) and in two Monte Carlo simulations. The solid simulation includes the 
s-wave amplitude 0.36 exp(i7r/4) (GeV)~^; while the dashed simulation neglects it. 
The known charm backgrounds are included in both cases. Only a small modulation 
is observed primarily in the tails due to the inclusion of the new amplitude. 
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